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Abstract 



It is demonstrated that any attempt to calculate the perturbative QCD contribution to the 
pion form factor requires the inclusion of intrinsic transverse momentum besides Sudakov form 
factors. For momentum transfers of the order of a few GeV the intrinsic transverse momentum 
leads to a substantial suppression of the perturbative QCD contribution. 
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There is general agreement that perturbative QCD in the framework of the hard-scattering 
picture is the correct description of exclusive reactions at asymptotically large momentum 
transfer. However, the applicability of this approach at experimentally accessible momentum 
transfers, typically a few GeV, was questioned [§, [3| . It was asserted that in the few GeV region 
the hard-scattering picture accumulates strong contributions from the soft end-point regions, 
rendering the perturbation calculation inconsistent. Recently, this statement was disproved by 
Sterman and collaborators [5], The essential point of their work is that the customarily 
neglected transverse momenta of the quarks as well as the Sudakov corrections are taken into 
account. Because of that the perturbative QCD contribution becomes self-consistent, even 
for momentum transfers as low as a few GeV, and seem to dominate form factors and cross- 
sections. For the electromagnetic form factor of the pion, the simplest exclusive quantity, it has 
been shown || that perturbative QCD can readily be applied for momentum transfers larger 
than 20Aqct>- It should be noted that the analyses of 0, [|, || involve no phenomenological 
parameter. 

However, it seems to us that an important element is missed in the analyses of 0, ||, namely 
the intrinsic transverse momentum dependence of the hadronic wave function. As we are 
going to demonstrate for the case of the pion form factor, the inclusion of that k±- dependence 
substantially suppresses the contribution from perturbative QCD. Although self-consistently 
calculated, this contribution then seems to be too small with respect to the data. We therefore 
suspect that there are sizable soft contributions to the pion form factor in the few GeV region, 
a conjecture already expressed in @, 0. 

In order to make our point of view clear it is useful to repeat the essential steps in the derivation 
of the hard-scattering formula for the pion electromagnetic form factor |]J . The starting point 
is the Drell-Yan formula in which the pion form factor is expressed as an overlap of the 
initial- and final-state light-cone wave functions 

Ftt(Q 2 ) = J *•(*!, k± + x 2 q) k ± ), (1) 

where Q 2 = q 2 is the momentum transfer, x\ is the usual longitudinal momentum fraction 
carried by the quark and k± its transverse momentum with respect to its parent hadron. The 
momentum of the antiquark is characterized by x 2 = 1 — x\ and — k±; \1/ is the wave function 
of the valence Fock state [} It satisfies the normalization condition 

f dx 1 d 2 k± -> . 2 

Contributions from higher Fock states are neglected in (|I]) since, at large momentum transfer, 
they are suppressed by inverse powers of Q 2 . As the inspection of ([!]) reveals, the large Q 
behaviour of the pion form factor is controlled by the tail of the wave function at large k±. The 
crucial point is that the tail of the wave function can be calculated, within perturbative QCD, 
from the soft part of the wave function by means of a homogeneous Bethe-Salpeter equation. 
Absorbing the perturbative kernel of the Bethe-Salpeter equation in a hard-scattering amplitude 
Tji, expression (]]]) can be converted into 

FM) = J / *5(lfc, h) T H {x x , y u Q, k ± , l ± ) k ± ). (3) 

Now \&o represents the soft part of the pion wave function, i.e. the full wave function with 
the perturbative tail removed from it. The perturbative tail of \1/ behaves as l/k\ for large k± 

4 Strictly speaking * is the S — 0, L z = component of the valence Fock state wave function. 



whereas \l/o vanishes as ~ 1/^1 or faster. To lowest order in perturbation theory, the hard- 
scattering amplitude is to be calculated from the one-gluon-exchange diagrams. Keeping the 
transverse momenta of the quarks, k± and l±, only in the gluon propagators, where they matter 
most 0, the hard-scattering amplitude reads 

rp i n Z T \ 16vr a s (/i)CV 

^, fc(? ,t ± ,/ ± ) = ___ ? _ E _, (4) 

where we have made use of the symmetry of the pion wave function under the replacement 
X\ <-> %2 (C invariance). In eq. (||) a s is the usual strong coupling constant to be evaluated at 
a renormalization scale [i and Cp (= 4/3) is the colour factor. 

At large Q one may neglect the k±- and /^-dependence in the gluon propagator as well; Th can 
then be pulled out of the transverse momentum integrals, and these integrations apply only to 
the wave functions. Defining distribution amplitudes (DA) by 

' f " <j>{xi,n) = f ^^y (xi,k±), (5) 



2^6 v ' J 16tt 
one arrives at the celebrated hard-scattering formula for the pion form factor 

F^ HSP (Q 2 ) = ^- J dxi dy x (j}*(y h /i) T H {xi,yi, Q, //) //), (6) 



which is valid for Q — > oo. 
The DA is defined such that 



o 



i 

dx\ 4>(xi, /i) = 1. (7) 



An immediate consequence of the definitions (|5]) and (0) is that the constraint from the 7r — > [iv 
decay is automatically satisfied: 

fn f dx\ d 2 k± -i 

^o[x 1 ,k±), (8) 



2VQ J 16tt 3 

where f w (= 133 MeV) is the usual n decay constant. The integral in ([5]) has to be cut off at 
a scale of order Q, which leads to a very mild dependence of the DA on the renormalization 
scale (QCD evolution). For most applications of perturbative QCD to exclusive reactions, the 
/i-dependence of the DAs can be ignored since one is only interested in a rather limited region 
ofQ. 

The neglect of the transverse momentum dependence in ([|) is a bad approximation in the 
end-point regions where one of the momentum fractions, Xi or yi, i = 1,2, tends to zero. 
This approximation is after all responsible for the theoretical inconsistencies mentioned in the 
introduction. Li and Sterman have suggested to retain the transverse momentum dependence 
of Th- Moreover, Sudakov corrections, suppressing the contributions from the dangerous soft 
regions, ought to be taken into account. In order to include these corrections it is advantageous 
to reexpress eq. ([3]) in terms of the Fourier transform variable b in the transverse configuration 
space 

F^ ert (Q 2 ) = ^^|r J^J 2 hi>l{yxXw)fH{x 1 ,y 1 ,Q^t)i> Q {x 1 -lw^ 

x exp [-S(x 1 ,y 1 ,Q,b,t)] (9) 



5 The effect of the /c^-dependence in the fermion propagators has been investigated by Li |6|; it leads to a 
reduction of the results for by about 10% for Q > 3 GeV. 

6 Since Tjj depends only on fcj^ + 1± there is merely one independent variable in the transverse configuration 
space. 



where the Fourier transform of a function / = f(k±) is defined by 



Kb) = J d 2 k ± e~ iU - f(k ± ). (10) 

The Fourier-transformed hard-scattering amplitude reads 

T H {x-i,y^Q,b,t) = 8a s (t) C F K (^/xTyiQ b), (11) 

where K is the modified Bessel function of order zero and t is the largest mass scale appearing 
in f H : 

t = Max{y/x^Q,w=l/b). (12) 

The factor exp [— S] in (Rl), termed the Sudakov factor, comprises the radiative corrections. The 
function S is given by [§J 

S(x u yi , Q, M) = E [s{xi, Q, b) + s( Vi , Q, b)} - |- In ^L^L , (i 3) 



The last term in (|13"D arises from the application of the renormalization group. The lengthy ex- 
pression for the Sudakov exponent s, which includes all leading and next-to-leading logarithms, 
is given explicitly in ||. The most important term in it is the double logarithm 

301 -/2Aqcd MIMqcd) 1 ' 

where £ is one of the fractions, x,i or y^, and (3\ = (33 — 2nf)/12. For small b, i.e. at small trans- 
verse separation of quark and antiquark, there is no suppression from the Sudakov factor]]. As b 
increases the Sudakov factor decreases, reaching zero at b = 1/Aqcd- F° r b larger than 1/Aqcd 
the Sudakov factor is set to zero. Owing to this cut-off the singularity of a s is avoided without 
introducing a phenomenological parameter (e.g. a gluon mass). For Q — > oo the Sudakov 
factor dumps any contribution except those from configurations with small quark- antiquark 
separation. In other words, the hard-scattering contributions dominate the pion form factor 
asymptotically. 

Li and Sterman have explored the improved hard-scattering formula (|9|) on the basis of cus- 
tomary wave functions, neglecting the QCD evolution and the intrinsic transverse momentum 
dependence of the wave functions. Their numerical studies have revealed that the modified 
perturbative approach is self-consistent for Q > 20Aq C d in the sense that less than, say, 50% 
of the result is generated by soft gluon exchange (a s > 0.7). In the few GeV region the values 
for F n as obtained by Li and Sterman are somewhat smaller than those provided by the hard- 
scattering formula @ and are, perhaps, smaller than the experimental values 0. The data 
may however suffer from large systematic errors p| . This prevents us from giving a definite 
conclusion about the agreement or disagreement between data and theory. 
The approach proposed by Sterman and collaborators 0, |§, certainly constitutes an enor- 
mous progress in our understanding of exclusive reactions at large momentum transfer. We 
believe, however, that in any practical application of that approach one has to allow for an in- 
trinsic transverse momentum dependence of the hadronic wave function, although, admittedly, 
this requires a new phenomenological element in the calculation. In order to demonstrate the 

7 Actually the Sudakov factor includes a small enhancement in the low-6 region. Li and Sterman [H neglect 
this enhancement and set any factor exp [— s(£, 6, Q)] to unity in this region and whenever £ < \/2/bQ. We 
apply this recipe too. 



importance of the intrinsic transverse momentum for the case of the pion form factor, we need 
realistic pion wave functions. 



The pion wave function: In accordance with (|5|), ([/]) and 
valence quark Fock state wave function as 



*o(zi,fci 



2^Q 



) we write the soft part of the 
xi) E(xi,fcj.), (15) 



the function E being normalized in such a way that 



16tt 3 



1. 



(16) 



The probability of the valence quark Fock state is given by (|2]) 0. A further constraint of the 
wave function comes from the charge radius of the pion. The (transverse) radius of the valence 
Fock state should be smaller than (or at best equal to) the full radius. Or, in terms of transverse 
momentum, we have a lower limit on the mean square transverse momentum that is defined by 

f^ 2 f dxi d 2 kj 



{kx') 



24 P qq J 167T 3 



ki\(f>(xx)\ 2 \^( Xl ,k ± )\\ 



(17) 



Actually, the root mean square transverse momentum (r.m.s.), (fcj_) 1//2 , should be larger than 
300 MeV f\ Lepage et al. ||, examining the process 7r° — > 77, derived yet another constraint 
on the pion wave function: 

U 2 r 

dxi4>(xi)Y>(xi, kx = 0). 



12 

For the kx- dependence of the wave function we assume a simple Gaussian 

E( Xl ,kx) = 167r 2 p 2 g( Xl ) exp (-g( Xl )(3 2 k 2 ± ) , (19) 

g(x\) being either 1 or l/x\X2. In the first case the wave function \l/o factorizes in x\ and 
kx, which manifestly breaks rotational invariance. This theoretical deficiency is unessential 
to the purpose of this article. The second, non-factorizing case is obtained from a harmonic 
oscillator rest frame wave function by equating the energy propagators in both the rest frame 
and the infinite momentum frame M. This particular fc^-dependence goes along with a factor 



exp 



-f3 2 m 2 / X1X2 in the DA. Here, m q is a constituent quark mass for which we choose a value 
of 330 MeV. The Gaussian fliTS| ) is consistent with the required large-fc^ behaviour of a soft 
wave function. For the DAs we try two versions, the asymptotic form ~ X1X2 and the CZ form 



X\X2(xi — X2) |Pl |12j. Thus all together we utilize four examples of pion wave functions: 



<Pa{x) 
4>b{x) 

<f) c ( X ) 



A a xix 2 

A b xiX2(xi - X2) 2 



9 
9 



A c xiX 2 exp 



21 



XiX 2 _ 

A d x 1 x 2 (xi - x 2 ) 2 exp 



1 



Pd 2 m 2 q 
X1X2 



9 



xix 2 
1 

xix 2 



(20a) 
(20b) 

(20c) 
(20d) 



8 The use of the soft part of the wave function in (g) instead of the full one, entails an insignificant 0(aj) 
error. 

9 Lepage et al. Q have proposed to determine the radius, and in turn the r.m.s. transverse momentum, 
from the derivative dF^/dQ 2 at Q 2 = which can be obtained through eq. (Q). Results obtained that way are 
consistent with our ones. 



Example b) emphasizes most the end-point regions, example c) least of all. In principle the 
DAs depend on the renormalization scale. Asymptotically, that is for Q — > oo, they all evolve 
into the asymptotic form (|20|a), which itself has no evolution. We are going to ignore the 
evolution in the main part of this article since it is a complication of minor importance. We 
will however comment below on the effects the evolution may cause. Finally, a remark is in 
order concerning the DAs (pO|c) and (f20|d). Such functions have been proposed by Chernyak 
and Zhitnitsky [|TT|] on the basis of the two lowest moments of the pion DA as derived from 



QCD sum rules. The CZ ansatz for the DAs is the subject of considerable controversy in which 
we do not want to enter 0. We merely consider (^c) and (pO|d) as examples of DAs strongly 
concentrated in the end-point regions. 

Our wave functions have one free parameter, the oscillator parameter /3, which we fix by requir- 
ing specific values for the r.m.s. transverse momentum, namely 350 and 250 MeV. In table 1 we 
compile the properties of our wave functions, namely the constants A, the oscillator parameters 
P, the probabilities of the qq Fock state P q q and the value of the right-hand side of eq. (fTH), 
the 7r — * 77 constraint. The value of 350 MeV for the r.m.s. transverse momentum turns out 
to be a reasonable value on all accounts. The corresponding radius is a little smaller than the 
measured charge radius of the pion and the constraint ([18]) is well satisfied. The probabilities 
of the qq Fock state also have plausible values. On the other hand, the value of 250 MeV for 
the r.m.s. transverse momentum is unrealistic. The radius is too large and ( |T8"D is strongly 
violated. We, however, keep this value for comparison. 

Numerical results: The Fourier transform of the fc^-dependent part of the wave function reads 

S( Il ,6) = te e X p(- i ^). (21) 

Li and Sterman Hj assume that the dominant 6-dependence of the integrand in eq. (|^) arises 
from the Sudakov factor and that the Gaussian in £ can consequently be replaced by 1. In 
order to examine that assumption, we compare in Fig. 1 the Gaussian with the Sudakov factor 
[see Cp~3|)l . We only display the case g = 1 since the other case, g = l/x±X2, is very similar. 
Obviously the intrinsic /cj_-dependence of the wave function cannot be ignored. For momentum 
transfers of the order of a few GeV the wave function damps the integrand in (|9]) more than 
the Sudakov factor. Only at very large values of Q does the Sudakov factor take over. What is 
more important at a given value of Q, the Sudakov factor or the wave function depends on the 
value chosen for Aq C d an d (k± 2 ) 1 ^ 2 . For the favoured values, 200 MeV for Aq C d and 350 MeV 
for the r.m.s. transverse momentum, the wave function plays the major role in the large b 
behaviour of the integrand in (|9]) for Q less than 20 GeV, whereas the Sudakov factor begins 
to dominate beyond that value, which is well above the experimentally accessible momentum 
transfer region. 

Numerical evaluations of the pion form factor through (H), using the wave functions given in 
(pQ"D, confirm the observations made in Fig. 1. For the two extreme DAs, the one that is most 
concentrated in the end-point regions, fl20| b), and the one that is least, (|20|c), we display the 
results in Fig. 2. The results for the other two examples, (|2"0|a) and (|2"0|d), lie in between 
the two extreme cases. As compared with the hard-scattering result, eq. the inclusion 
of the Sudakov factor leads to some suppressions as already noticed in ||. The intrinsic k±- 

10 For instance, Braun and Filianov also employing QCD sum rule techniques, found a value of 1.2 ± 0.3 
for the pion DA at x = 1/2, in clear contradiction with the DAs of the CZ type. It has been pointed out by 



Radyushkin [|| that the finite size of the vacuum fluctuations had not been considered in |11 . Taking these 
finite-size effects into account, one finds values for the moments of the pion DA that are much closer to those 
obtained from the asymptotic DA than the values quoted in O] . 



dependence of the wave functions provides additional suppression which amounts to about 
50%, in the case fl20"|b), at Q — 2 GeV and is still substantial at 7GeV. As expected the total 
suppression obtained from both the Sudakov factor and the intrinsic fcj_-dependence is stronger 
for the wave function (pO|b) than for (|20|c). One may regard the Sudakov factor as a part of 
the wave function, in fact as its large-6 behaviour. In this case one would not consider the 
product of the two functions, exp [-S] and S. Rather, to an admittedly crude approximation, 
one would take either the Gaussian in £ or the Sudakov factor, whichever is the smaller at a 
given value of b. Fortunately this recipe yields only minor differences with the first case. 
The strong suppression due to the intrinsic fc^-dependence of the wave function reveals that the 
approximation made in the derivation of the hard-scattering formula (|^), namely the neglect of 
the /c^-dependence in the hard-scattering amplitude (f|), is invalid at small momentum transfer. 
The reason for this can easily be understood. Using, for instance, a wave function of the type 
(|15],|19D with g = 1 and, for the sake of the argument, replacing a s in @ by an average value, 
one can convert (131) into 



Q 2 F W (Q 2 ) = ^U 2 a s C F \ U dze z E x (z) f %{y/u)<t>{z/y) 
6 Jo Jz y 



(22) 



where u = Q 2 (5 2 /2 and E\ is the exponential integral. The dominant contribution comes 
from the upper limit of the z integration. Replacing E\ by the leading term of its asymptotic 
expansion 

ze z E 1 (z)^l-- + l T ... : (23) 

z z A 

one obtains the hard-scattering contribution (^|). But the higher-order terms provide substantial 
corrections to it. The wave function (|20| a), for instance, leads to 



Q 2 FAQ 2 ) = Q 2 f„ hsp {Q 



-?+■■■' 



(24) 



In Fig. 2 we have also shown the available data for the pion form factor at large momentum 
transfer PI. This is not meant as a serious comparison between data and theory since, as we 



mentioned above, the data may suffer from large systematic errors [ 1Q| . Nevertheless, it seems 
to us that the perturbative contribution (^|), including the Sudakov corrections and the intrinsic 
fc_i_-dependence of the wave function, is too small with respect to the data, perhaps by a factor 
of 2 or so. One may suspect the QCD evolution of the DAs, ignored by us up to now, to be 
responsible for that discrepancy. Let us estimate whether this explanation holds or not. The 
evolution equation for the pion DA has been derived in [|J]. For our example (|20|b) it reads 



(f) b (xi,w) = a (Xi) 



l + -C 3 2 /2 (x 1 ~x 2 ) 





50/81" 


a s (w 2 ) 


.a s (^o). 





(25) 



where the C^ 2 (C) are the Gegenbauer polynomials, which are the eigensolutions of the evolution 
equation. According to || the evolution parameter w is taken to be 1/6. Inserting (^) into 
eq. (H) and using fi = 0.5 GeV we re-evaluate the pion form factor and find that the 
evolution reduces even further the predictions for F n (by about 30% for Q between 3 and 4 
GeV). Similar results are obtained for the other examples of wave functions, (pO|c) and (|20|d). 
This is in agreement with the findings reported in |i~2"|. 



Thus we are tempted to suppose that higher-twist contributions, besides the /cj_-dependence, 
are quite strong in the few GeV region. Sources of such higher twists are contributions from 
higher Fock states, L / components in the valence Fock state wave function or genuine soft 



contributions. 

An example of genuine soft contributions is provided by the starting point of our considerations, 
namely the Drell-Yan formula ([I]). Up to now we have only considered those contributions from 
(P which arise from the perturbative tail of the wave functions. The size of the ignored soft 
contributions until now can easily be estimated for our wave functions (EDI) . One finds 



2 

F, soft = y f 2 P 2 J dxt g{x 1 ) <p 2 { Xl ) exp (-g( Xl ) (3 2 x 2 2 Q 2 /2) . (26) 

The integral is dominated by the region near x\ = 1. Other x± regions are strongly damped by 
the Gaussian. Hence F n sensitively reacts to the behaviour of the DA for x\ — > 1. Numerical 
evaluations of (|26|) for our wave functions (|2~0"D show that the examples fl20la) and (|20|c) provide 
soft contributions of the right magnitude to fill in the gap between the perturbative contribu- 
tion (|) and the experimental data Q (Q 2 F* oft ~ 0.25 GeV 2 at Q = 2GeV). As required by 
the consistency of the entire picture, F^°^ decreases faster with increasing Q than the pertur- 
bative contribution. The exponential exp(— (3 2 m 2 / xix 2 ) in the example (|20|c) turns out to be 
quite important since it is effective in the end-point regions. It reduces the size of Q 2 F^ ^ 
substantially and alters its asymptotic behaviour from a 1/Q 2 behaviour to an exponentially 
damped decrease; Q 2 F^ ^ becomes equal to the perturbative contribution at Q ~ 5 GeV. The 
soft contributions are also subject to Sudakov corrections. An estimate of these corrections on 
the basis of eq. (|T3p reveals, however, that they only amount to a few per cent for the examples 
(pO|a) and (pO|c). The wave functions (pOfc) and (pQ|d) provide very large soft contributions be- 
cause of their strong concentration in the end-point regions. Therefore, these examples appear 
to be unrealistic [see also the discussion after eq. (|20|)1. Soft contributions of the type ( [26] ) 
have extensively been discussed in 0. Our results for F* ^ are similar in trend, but smaller 
in size than those presented in [[|. There are three reasons for the differences in the size of 
the soft contributions: Isgur et al. @ assume that P q q = 1, they neglect the exponential 
exp(— f3 2 m 2 / x\x 2 ) and utilize DAs that are proportional to ^X\X 2 instead of ~ X\X 2 . Strong 
soft contributions to the pion form factor have also been obtained with QCD sum rules [[3]. 
Finally we mention L ^ components of the wave function as another source of higher-twist 
contributions. For instance, the valence Fock state component of the pion may be expressed by 

/ ^6^7! U + m) [*°( a?1 ^ J -) + ¥±^iM)] 75, (27) 



where p denotes the pion's momentum and m its mass. Ralston and Pire [HJ consider an 
even more general ansatz for the valence Fock state wave function of the pion. Thus one may 
have additional fc^-dependent contributions to the pion form factor, which are not necessarily 
small. Perhaps one obtains a better description of the pion form factor with such contributions 
although at the expense of the introduction of a new, a priori unknown, phenomenological func- 
tion. Contributions from higher Fock states suffer from the same disadvantage; also in this case 
new phenomenological functions have to be introduced. Such extensions of the hard-scattering 
model are certainly not very attractive but are, perhaps, necessary. 

Summary: On the basis of our numerical studies we conclude that the intrinsic k±- dependence 
of the wave function has to be taken into account for a reliable quantitative estimate of the 



n The Drell-Yan formula (Q) has been obtained by a transverse boost to the infinite-momentum frame. Since 
we know the hadronic wave function only approximately, it may happen that other boost directions provide 
different results for the pion form factor. For a discussion of such eventual violations of rotational invariance, 
we refer to §, Q. 



perturbative QCD contribution to the pion form factor. We are aware that this introduces 
a new phenomenological element into the calculation. The disadvantage is, at least partially, 
compensated by the fact that the inclusion of the intrinsic fc^-dependence renders the per- 
turbative contribution even more self-consistent than the Sudakov suppression already does. 
Applying the criterion of self-consistency as suggested by Li and Sterman (see above), we 
can conclude that perturbative QCD begins to be self-consistent for Q between 1 and 2 GeV 
(for (k 2 ^) 1 / 2 = 350 MeV). The value for Q at which self-consistency sets in depends on the wave 
function. It is larger for the end-point concentrated wave functions (|20"|p,d) than for the other 



examples (20a,c). However, the perturbative contribution, although self-consistent, is presum- 
ably too small with respect to the data. It thus seems that other contributions (higher twists) 
also play an important role in the few GeV region. We see no argument why our observations 
should not apply as well to other exclusive observables, such as the nucleon form factor. 
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Table 1: Properties of the pion wave functions (20). 



DA 
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(3 2 [GeV- 2 } 


p. 


R.h.s. of ( 
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(3 2 [GeV- 2 } 


p. 


R.h.s. of (PD 


a) 


6 


4.082 


0.284 


0.947 


6 


8.0 


0.557 


1.856 


b) 


30 


4.082 


0.338 


0.947 


30 


8.0 
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10.07 
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1.044 
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1.787 
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2.006 


d) 


50.05 


0.571 


0.357 


0.983 


83.59 


1.247 


0.721 


1.915 








= 350 MeV 






= 250 MeV 



Figure captions 



Fig. 1: The Sudakov factor, evaluated at x\ = yi = 1/2, and the Gaussian exp ( — 6 2 /2/3 2 ) [see 



eq. (21f)1 as functions of the transverse separation b. The Gaussian is shown for a r.m.s. 
transverse momentum of 350 MeV (short-dashed line) and 250 MeV (long-dashed line). 
The Sudakov factor is evaluated for Aq C d = 200 MeV (solid lines) and 100 MeV (dotted 
lines) for Q 1 = 2 GeV, Q 2 = 5 GeV and Q 3 = 20 GeV (a* = x 2 = 1/2). 

Fig. 2a: The pion form factor as a function of Q 2 evaluated with the wave function (pO|b) and 
A.qcd — 200 MeV. The dash-dotted line is obtained from the hard-scattering formula with 
an a s cut-off at 0.5 and the dashed line from (^|) ignoring the intrinsic fc^-dependence. 
The solid and the dotted lines represent the complete result obtained from @ taking into 
account both the Sudakov factor and the intrinsic fcj_-dependence for (k^) = 350 MeV 
and 250 MeV, respectively. Data are taken from [[| (o 1976, • 1978). 



Fig. 



2b: As Fig. 2a but using the wave function (pO|c). Note the modified scale of the abscissa. 
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